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It is found that Grover's quantum search algorithm is not robust against phase inversion and 
Hadmard transformation inaccuracies. Imperfect phase inversions and Hadmard-Walsh transforma- 
tions in Grover's quantum search algorithm lead to reductions in the maximum probability of the 
marked state and affect the efficiency of the algorithm, even in the absence of decoherence. Given 
the degrees of inaccuracies, we find that to guarantee half rate of success, the size of the database 
should be in the order of 0{-h), where 5 is the uncertainty. 
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O ■ I. INTRODUCTION 

o 

\^ ' Grover's quantum search algorithm is a remarkable achievement in quantum computing |l[]. There have been 

intensive interests in Grover's quantum search algorithm recently H-K] It uses only two simple gate operations, 
the controlled phase rotations and Hadamard transformations. It has been successfully demonstrated in solution 
NMR bulk quantum computers with a few qubits |0,pi. However, the inevitable quantum state decoherence and gate 
inaccuracies can introduce errors ||18| , |l9| , which accumulate through the computation and make long computation 
unreliable. While, in order to find out the marked state with high probability, it still requires exponential number of 
f~^ . iterations. Then, the error probability of the complete algorithm may be as exponentially large as the error probability 
C^ ' of each iteration. In other words, even with small imperfection per step, large scale quantum search may be difficult. 
T-H , Fortunately, recent study in quantum error correction shows that in principle, whenever the noise rates are below a 

0^ '■ constant threshold, an arbitrary long quantum operations can be performed reliably through fault-tolerant quantum 
^^ computation [ pO| . Experimentally, different types of faults can occur with different rates and will affect the efficiencies 

C^ of the algorithm differently. For example, the effect of quantum state decoherence and operational errors on the 
Oh! efficiency of quantum algorithms have been studies in pl[ with ion trap quantum computers. A good understanding 
4_j ■ of the effect from different noises on the algorithms can help us look for specific potential physical realizations of 
^ ^ quantum computers. 

In this paper, we address the problem of influences of imperfect gate operations in the quantum search algorithm, 
in the absence of decoherence and error corrections. We will show that systematic phase mismatching and random 
errors in the Walsh-Hadmard transformation lead to exponential reduction in the maximum success probability when 
n is linearly increased. Therefore, if we can not avoid them completely, to ensure a large success rate in a quantum 
searching machine, the size of the database should be limited. This limitation is due to the intrinsic vulnerability 
5_^ ■ of the algorithm to imperfect gate operations. In designing a quantum searching machine, this limitation should be 
Ci ' taken into account. 

The paper is organized as follows: Section 2 is devoted to the description of different error models in phase mis- 
matching and the corresponding simulation results. In section 3, we present the consequence of imperfect Hadamard 
transformation. Section 4 gives a short summary. 

II. EFFECTS OF IMPERFECT PHASE INVERSIONS 

Grover's algorithm consists of essentially four steps in an iteration pj: (1) a Walsh-Hadamard transformation 
U ~ W; (2) a phase inversion of the prepared state I7), usually I7) = |0), I^ — I — 2|7)(7|; (3) a phase inversion of 
the marked state \t) , Ir = I — 2|t)(t| ; and (4) an inverse of the Walsh-Hadamard transformation U^^ = W {W is 
self-inverse.). The operator for one Grover iteration is Q = —I~^U^^IrU. 
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In this section, we focus on the imperfection in phase inversions and therefore choose U to be the ideal Hadamard 
transformation. We consider the imperfections in the phase inversion to be systematic, such that 

I,=I-{l~e^')\l){ll 
/,=/-(l-e*'^)|r)(r|, 

(1) 

where 9 — -k + 9q,lp — tt + lpq with ^o and tpo constant and small. When 9q — ^po ^ 0, we recover the original Grover's 
algorithm. The generalized quantum search algorithm is a rotation in a 2-dimensional space spanned by I7) and \t). 
In the following two orthonormal basis 



Il) = 



(|7)-C/.,C/-») 



|2) = C/-i|r). (2) 



with Ur-y = {t\U\^) = l/\/N, the operator Q is represented by 






21 ) (3) 



Let 5 = 9 — ip = 9q — LpQ. It has been shown that to construct an efficient quantum search algorithm, 9 and Lp must 
equal to one another |2^j2^. However due to imperfections in gate operations, this phase matching requirement can 
not be strictly satisfied. In the following, we show that nonzero constant 5 results in exponential reduction in the 
maximum success probability of Grover's algorithm asymptotically. 

Since both ^0 and ^pq are small, dropping off an overall phase, we approximate Q as 

Q = cos (5/ + i sin (5cr^ + 1(5' dy + o(/3') , (4) 

where ax, o-y and cr^ are Pauli operators and / is the identity operator in dimension 2. /?' == 2/3 + O{9ol3) —2(3 + o{/3) 
with (3 = ^ ~ . For small S , we can further simplify operator Q as, 

with G = sinSa^ + P'uy. Using G^ == ((5^ + /3'2)/, we obtain 



cosjA + z^^ ^^^ 



(5) 



with A ^ yfW+W^. Then, starting from the prepared state I7) = ^\ - |C/^^|2|1) +[/^^|2) = cos/3|l) +sin/3|2) w |1), 
after j number of iterations, the norm of the amplitude of the marked state in the quantum computer is 



I. 
A 



|B,|«^sin(jA). (6) 



and the maximum probability of the marked state in the algorithm is 

/3'2 



<1- (7) 



Therefore, for large N , Grover's algorithm is efficient only when (5 = 0. When (5 7^ 0, we find 



(8) 



decreases linearly with N or exponentially with n = log2 N. This concludes our proof that systematic 
phase mismatching results in exponential reduction in the success probability and consequently gives an upper bound 
on the size of the database. If half rate of success is required, that is Pmax > 1/2, N cannot exceed 8/S'^. 

So far, we have assumed that the errors in the phase inversions are systematic such that S is constant. We now 
extend this simple model (EMI) to another two error models. The second error model (EM2) assumes S in each 



step is a Gaussian random variable with mean So = and standard deviation s. Such an error is conventionally 
defined as random error. Finally, we let 6 he a Gaussian random variable with mean So j^ and standard deviation s 
(EMS). The exact effect of EM2 and EM3 are difficult to compute analytically due to their randomness. Hence, we 
only present the simulation results. We vary n = Log2N and run the algorithm with sufficient number of iterations 
so that a maximum probability is found. Since S in EM2 and EMS are random variables, we adopt the random 
sampling techniques in the simulation. The relationships between the maximum success probabilities and the size of 
the database are shown in Fig. 2 and Fig. S for EM2 and EMS respectively. For comparisons, we also provide the 
simulation result from EMI in Fig.l. 

Our simulation results are consistent with mathematical predictions. First, both systematic and random errors 
cause reduction in the maximum probability. Second, the success probability drops quickly after a transition point, 
which is determined by the error parameter Sq and s. When n is large, the probabilities decreases exponentially. Third, 
the different effects of systematic errors and random errors also meet our expectations. Mathematically, systematic 
errors cause the error amplitudes to grow exponentially with the number of gates applied; while the random errors 
cause the error probabilities to grow linearly This difference is clearly demonstrated in our simulation results. Fig. 2 
shows that random errors give a much larger transition point than systematic errors. Fig. S shows that the average 
success probability from EMS is nearly identical to EMf except some small fluctuations. 

It is shown in this section that systematic errors in the phase inversions lead to reduction in the maximum probability 
of finding the marked state. Random erros also affect this successfl rate, but in a lesser degree. In practice, we should 
make Sq as small as possible. However, due to imperfection, nonzero Sq occur inevitably, due to imperfection, nonzero 
d occur inevitably. For instance, systematic errors arises from imperfect calibration and inhomogeneity in the radio 
frequency pulses in NMR realization. Random errors are always present in a realistic environment. These errors will 
reduce the maximum probability of the algorithm. To make an estimate on the combined effect of systematic and 
random errors(EMS), we assume that random errors affect the algorithm just like the systematic errors. Then we can 
treat A = 25 as the uncertainty due to both systematic errors and random errors and use this to derive an upper 
bound for the size of a quantum database: any phase inversion operation is imperfect, there is an uncertainty, and this 
uncertainty sets an upper bound on the size of the database N . For half rate success, the dimension of the database 
should be less than ^. 

III. IMPERFECT HADAMARD TRANSFORMATION 

Hadmard- Walsh transformations are also subject to errors. To study the effect of the imperfect Hadmard- Walsh 
transformation, let's take 5 = in eqn. (ra). Then the maximum probability for finding the marked state is approxi- 
mately sin^(2j7?) for perfect unitary transformation. For perfect Hadamard- Walsh transformation, /3 — arcsin(|[/T-^|), 
and \Ur-y\ = y^l/N. For systematic errors in the Hadmard- Walsh transformation, the matrix elements of U is no 
longer equal to y/l/N. If lUr-y] is larger than \/ jt, then the algorithm will require less steps in reaching the desired 



state compared with the standard Grover's algorithm. If it is smaller than \/jt, the algorithm will require more 
steps of iteration. In this case, the searching algorithm still can give a probability quite close to unity. But if one 
makes a measurement at the normal optimal number of iteration, one will gets a reduction in the successful rate. 
This difficulty can be overcome by using the algorithm several times with measurements made around the optimal 
iteration which is similar to the one used in Ref. |lj] . 

Here, we can give a simple interpretation why Grover's algorithm is optimal. The rigorous proof has been given 
in Refs. ||ll|]. Grover's algorithm can be seen as a rotation of the state vector in a 2-dimensional space span by \t) 
and 7). Each iteration rotates an angle A = /3' = 2 sin(6'/2)/3. 6 = (p = tt gives the largest angle 2/3 — 2arcsin(|[/T--y|. 
So one has to choose phase inversions. As for the unitary transformation U, at first glance one maybe attempted to 
think that a larger {UtjI will constitute a faster search algorithm. However, since U is unitary, its matrix elements 
satisfy the normalization relation J^t |f^T7p — 1, where r runs through all the N basis states. The mean value of 

the matrix element is \/ jt- If some of the matrix elements are larger than the average, some other matrix elements 
will be less than this average. In other words, while making the search for some marked states in less steps, the 
modified algorithm has to search the rest of the basis states in more steps. In contrast, the original Grover algorithm 
searches all possible marked state with the same optimal number of iterations. Together with its simpleness and easy 
implementation, the Walsh-Hadamard transformation lend itself the best choice. 

We discuss the effects of random errors in the Walsh-Hadamard transformation in a simple model. In this case, the 
algorithm is no longer a simple rotation in 2 dimensions. Though in each iteration, the operator can be approximately 
written as 



/ cos/3 sin^\ 

y~sinp cosp y ^ ^ 

the basis states in each iteration has been changed, that is, the 2 dimensional space in each iteration is no longer 
the same. This is apparent by inspecting the expressions in cqn. (^. Suppose in the first iteration, the unitary 
transformation is U and in the following iteration, the operator becomes V. Then after the first iteration, the state 
vector of the quantum computer is 

IV'i) = cos/3|l) - sin/3|2) « cos/3|l') - sin^C/-V|2'), (10) 

where |2') = ^"^|r). Because C/ ^ V, [/-ly is no longer the identity operator. Expanding J7-iy|2') = ((7-^)22 12') + 
..., we see that the Grover search operator acts only on the subspace span by |1') and |2'), and the other terms are 
leaked out the 2 dimensional space. To make an estimate, let's assume that in each iteration, (t/~^y)22|2') w (1 — (5i)|2) 
+ higher order terms. Then in this model, the matrix for a Grover search operator becomes 

_f cos/3 sin/3(l -(5i) \ , . 

-sin/3 cos/3(l-(5i) ) ' ^ ^ 



Starting from initial state I7) w |1), after j iterations, the amplitude of the state |2) becomes 

(12) 



(l-:^5i)sin(j/3) 



where only first order in ^1 is retained. With optimal number of iterations, j « ^ /^ , sin(j/3) « 1, the successful rate 

is 



P. 1-^ .1-!^. (13) 

For half success rate, one must have N < —$ri, which is similar to the limitation on the size of the database in the 
phase inversion inaccuracies. However, the mechanism is different. Here the random errors play a more important 
role than the systematic errors, whereas in the phase inversion case, it is just the opposite. 

IV. SUMMARY 

In summary, we find that the dominating gate imperfection in Grover's algorithm is the systematic phase mis- 
matching and the random errors in the Walsh-Hadmard transformation. Using the results obtained in this work, it is 
easy to understand the simulating results of Ref. pl[. In Fig. la of pl| is the results with only random errors(both in 
phase inversions and Hadmard transformation), we see that the peak in the probability curve drops down as random 
errors grow. But the position of the peak is relatively fixed. Random errors in the phase inversions does not affect 
the algorithm very seriously. Random errors in the Hadmard transformation reduce the maximum probability. The 
optimal iteration number remain more or less the same. When there is only systematic errors as shown in Fig. 2b of 
Ref. [ pl| , we see a drop in the maximum probablity and also a shifting of the peak position to the left, as is shown in 
Fig. 2b of Ref. pTI. The drop in maximum probability is caused by phase mismatching. The shift of the peak position 
is due to the systematic errors in the Walsh-Hadmard transformation. 

These gate inaccuracies set an upper bound on the size of the database. We estimate that the upper bound 
is inversely proportional to the quadrature of the uncertainty in the phase mismatching or in the Walsh-Hadmard 
transformation. In real quantum computation, imperfect gate operations exist all the time at constant rate while 
decoherence increases rapidly with computing time. At the early stage of a quantum computation, gate imperfection is 
dominant in affecting a quantum algorithm. As the computation continues, decoherence increases and then dominates. 
Suitable quantum correction codes and in particular fault-tolerant quantum computation can reduce the decoherence, 
and ease the stringent requirement on gate accuracies. The limitations on the quantum datasize can then be greatly 
relieved. 
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Figure captions Fig.l. EMI with So = 10"^^ 10"^^ 10"''. 
Fig.2. EM2 with Sq = 0,s = lO^^ and A = /?' 
Fig.3. EMS with So = 10"^ lO^^, lO"*, s = lO^^ 
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